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Abstract. The sibling species Littorina plena and L. scutulata show considerable intraspecific varia¬ 
tion in shell morphology, with wide overlap in shell characters between species. Measures of single 
characters (length-to-width ratio, apical angle, number of whorls, degree of tessellation, and presence 
of one spiral amber band in the shell aperture) are not sufficient for accurate separation and classifi¬ 
cation of specimens. However, discriminant functions that use statistically derived combinations of 
several easily measured shell characters allow approximately 92% of the snails from the exposed coast 
near Bodega Bay, California, to be correctly classified, and significantly improve morphological dis¬ 
crimination of these sibling species. 

The shell morphology of Littorina scutulata and L. plena varies geographically and between habitats. 
Variability in morphological characters may be as great within a species as between species when 
individuals are considered from several sites. Only 47% of individual L. plena are correctly classified 
when the discriminant functions derived for open coast populations are applied to snails from nearby 
protected habitats. Although discriminant analysis is a powerful technique for distinguishing between 
morphologically similar species, such variation in morphology emphasizes the need to derive and eval¬ 
uate discriminant functions for specific populations. 

INTRODUCTION 

Taxonomic status has been the subject of study for sev¬ 
eral species of marine littoral gastropods in the genus Lit¬ 
torina. Intraspecific variation in shell characters and wide 
overlap between species have at times made morphological 
identification of species complicated and ambiguous, and 
taxonomists have had to rely upon other sources of infor¬ 
mation to detect or confirm the presence of sibling (mor¬ 
phologically similar) species of Littorina. Differences in 
mode of reproduction (Heller, 1975a), morphology of 
genitalia (Whipple, 1965; Heller, 1975a; Goodwin & 

Fish, 1977), characteristics of spawn (Whipple, 1965; 

Borkowski & Borkowski, 1969), and allele frequencies 
at certain gene-enzyme loci (Heller, 1975a; Ward & 

Warwick, 1980; Wilkens & O’Regan, 1980; Mastro 
et al. , 1982) have been used to establish the existence of 
separate species in cases where morphological data on shell 
characters alone were insufficient. 

Within populations, many Littorina species exhibit con¬ 


siderable variation in shell size, shape, sculpture, color, 
and color pattern (for examples: Borkowski, 1975; Raf- 
faelli, 1979). Further morphological variation within a 
species can be attributed to environmental gradients in 
temperature (Hughes, 1979), substratum (Heller, 
1975b), and exposure to wave action (Struhsaker, 1968; 
Newkirk & Doyle, 1975; Hylleberg & Christensen, 
1977). Distinguishing between sympatric sibling species 
on the basis of shell morphology alone is tenuous when 
such intraspecific variability is present. 

Where sympatric sibling species do occur, morpholog¬ 
ical classification of individuals is improved by removing 
interhabitat variation in shell characters. Separate diag¬ 
nostic characters can be selected for each habitat in which 
the species coexist. If necessary, statistically derived com¬ 
binations of shell characters can be employed to discrim¬ 
inate between sibling species in the same habitat. 

Littorina plena Gould, 1849, and Littorina scutulata 
Gould, 1849, occur in sympatry along the west coast of 
North America. They are abundant microalgal grazers on 
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500 m 

Figure 1 

Location of study area on the biological reserve of the University 
of California Bodega Marine Laboratory: triangles indicate open 
coast sites and circles indicate protected harbor sites from which 
Littorina specimens were collected. 

rocks and pilings in sheltered bays as well as on rocky 
shores of the exposed open coast. Reproductive distinc¬ 
tions detailed by Murray (1979) and patterns of genetic 
differentiation described by Mastro et al. (1982) have 
verified the taxonomic status of L. scutulata and L. plena 
as separate sibling species. Both species, however, possess 
variable shell characteristics, and consistent differences in 
shell morphology have yet to be defined completely for 
these snails. 

Murray (1982) uses the statistical techniques of dis¬ 
criminant analysis for morphological classification of Lit¬ 
torina plena and L. scutulata from a few locations along 
the west coast of the United States. However, satisfactory 
species characterization of the morphology of Littorina is 
likely to require control of local environmental (habitat) 
variables. This study describes in detail the morphology 
of snails from several populations near Bodega Bay, Cal¬ 
ifornia, using techniques similar to those of Murray 


(1982), and emphasizes the potential hazards of using 
discriminant functions for organisms that show wide¬ 
spread geographical and habitat variation in morphology. 

MATERIALS and METHODS 
Exposed Coast Populations of Littorina 

Specimens of Littorina plena and L. scutulata were col¬ 
lected from exposed rocky shores northwest of Horseshoe 
Cove on the biological reserve of the University of Cali¬ 
fornia, Bodega Bay, California (Figure 1). Individuals 
were separated by species and sex according to the char¬ 
acters described by Murray (1979) and Mastro et al. 
(1982). Male snails were identified by penis morphology 
and female snails by the type of egg capsule produced 
while isolated in individual containers. Only identifiable 
snails (spawning females and males with shell lengths 
greater than 4 mm) were used in this study of shell mor¬ 
phology. 

Statistical analyses of morphological data were con¬ 
ducted on several shell characters, which included shell 
length, shell width perpendicular to the long axis, apical 
angle, number of whorls, degree of tessellated color pat¬ 
tern, and presence of one amber spiral band in the aper¬ 
ture. These characters were chosen for their ease of de¬ 
termination and for their initial promise in distinguishing 
Littorina plena from L. scutulata. Linear dimensions were 
measured to the nearest 0.01 mm with dial calipers held 
under a dissecting microscope and the apical angle was 
determined to the nearest degree with a protractor. The 
degree of tessellation and the appearance of the amber 
band were scored as either (0) absent, (1) obscure, or (2) 
distinct. All measurements were made on live animals 
with the aperture oriented upwards (Figure 2). 

Discriminant Analyses for Exposed 
Coast Populations 

An optimal set of shell characters for classifying ex¬ 
posed coast specimens of Littorina plena and L. scutulata 
was obtained by discriminant analysis. Discriminant 
analysis is a procedure that weights and linearly combines 
characteristics for which groups are expected to differ; the 
result is a discriminant function that maximizes statistical 
separation of the groups. It is then possible to identify 
those characteristics that contribute most to differentiation 
between the groups, and to derive functions that classify 
cases of unknown membership. 

Discriminant analyses were performed on several sub¬ 
sets of the data, using a stepwise method of entering and 
removing shell characters as discriminating variables. The 
computer program for discriminant analysis obtained from 
the Statistical Package for Social Sciences (NlE et al., 1975) 
began by selecting that character that best distinguished 
Littorina plena from L. scutulata. A second character was 
then selected as the variable best able to improve the pow¬ 
er of discrimination when combined with the first char- 
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Figure 2 

Morphological measurements made on each shell. L = shell 
length; W = shell width; A = apical angle. (Shell shown with 
amber band present inside the aperture.) 




acter. The third and subsequent characters were selected 
similarly in order of their explanatory power until all 
characters were selected. At each step, characters could be 
removed if they acted to reduce discrimination in combi¬ 
nation with more recently selected characters. Standard¬ 
ized (weighting) coefficients were eventually produced that 
indicated the relative importance of each of the characters 
finally chosen. Cooley & Lohnes (1971) and Tatsuoka 
(1971) have described the mathematical techniques in¬ 
volved. 

Snails from exposed coast populations were used to de¬ 
rive discriminant functions for Littonna plena and L. scu¬ 
tulata. From samples collected monthly during the period 
September 1980 to December 1981, a subsample of 275 
L. plena and 253 L. scutulata was selected randomly for 
discriminant analyses. Smaller sample sizes were used for 
some analyses, and individuals selected were restricted to 
shell lengths less than 13 mm because all snails larger 
than 13 mm were L. scutulata. 
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LENGTH/WIDTH RATIO 


Figure 3 

Frequency histograms of length-to-width ratios for Littorina plena 
(n = 275) and L. scutulata (n — 253). Distributions for the two 
species are not statistically different. Kolmogorov-Smirnov two- 
sample test: DM AX = 0.1001; P > 0.10. 


and 33 L. scutulata were obtained from the exposed coast 
habitat to examine the usefulness of the discriminant func¬ 
tions over long periods of time. The capacity of the dis¬ 
criminant functions to classify individuals from a different 
habitat was measured on a sample of 57 L. plena and 22 
L. scutulata obtained from the protected harbor habitats. 
Harbor snails were identified and measured using the same 
methods as described for exposed coast snails. 


Suitability of Derived Discriminant Functions 

Specimens were collected between May 1984 and Oc¬ 
tober 1984 from the exposed coast study site and from 
additional sites inside a protected harbor located 1.5-2.0 
km from the exposed coast site (Figure 1). These samples 
were used to test the general applicability of the previously 
derived discriminant functions. Thirty-four Littorina plena 


RESULTS 

Exposed Coast Populations of Littorina 

Littorina plena and L. scutulata from the exposed coast 
exhibited considerable overlap in all shell characters ex¬ 
amined in this investigation. Frequency histograms of 
length-to-width ratios (Figure 3) failed to reveal any sta- 
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PRESENCE OF AMBER BAND PRESENCE OF TESSELLATIONS 

Figure 4 

Frequency histograms of shell characters for Littorina plena and L. scutulata. n p = sample size for L. plena ; n s = 
sample size for L. scutulata. Species differences in all shell characters are significantly different (Kolmogorov- 
Smirnov two-sample tests; all P < 0.01). A, apical angles (rc p = 216; n $ = 174; DMAX = 0.2208); B, number of 
whorls (rc p = 333; n s — 315; DMAX = 0.4783); C, occurrence of amber band in the aperture (0 = absent; 1 = 
obscure; 2 = distinct) (n p — 280; n s = 272; DMAX = 0.6692); D, occurrence of tessellated color patterns on shell 
(0 = absent, 1 = obscure; 2 = distinct) (n p = 92; n s = 122; DMAX = 0.4073). 


tistical differences in general shell shape between these 
two sibling species. However, specimens of L. scutulata 
attained a larger size than those of L. plena (snails with 
shell lengths greater than 13 mm were invariably L. scu¬ 
tulata). Moreover, shells from these populations of L. plena 


and L. scutulata displayed statistically significant differ¬ 
ences in the apical angle, number of whorls, degree of 
tessellation, and presence of an amber band in the aper¬ 
ture (Figure 4). 

Individual Littorina plena generally had larger apical 
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angles, and the population as a whole showed greater 
variation (F = 1.95; d.f. = 215, 173; P < 0.01), than L. 
scutulata. The average apical angle (mean ± SD) for L. 
plena was 62.6 ± 6.0 while apical angles averaged 60.2 ± 
4.3 for L. scutulata. Figure 4A presents sample sizes and 
distributions of apical angles for each of the sibling species. 

Shells of Littorina plena usually had fewer whorls than 
those of L. scutulata (Figure 4B). Only 20% of the L. plena 
specimens examined (n = 335) had four or more whorls, 
while 68% of the L. scutulata shells ( n =315) had at least 
four whorls. 

Littorina plena shells were more likely to possess an 
amber band in the aperture (Figure 4C), but less likely 
to show tessellations (Figure 4D) than those of L. scutu¬ 
lata. In 87% of the L. plena specimens (n = 280) the band 
was present to some extent; the band was observed in only 
20% of the L. scutulata shells (n — 272). Thirty-seven per¬ 
cent of the L. plena shells inspected (n — 92) displayed at 
least some degree of tessellation; in contrast, a tessellated 
pattern was present on 76% of the L. scutulata shells (n = 
122), and this pattern was more often very prominent than 
on L. plena shells. 

Discriminant Analyses for Exposed Coast Populations 

Although the differences between Littorina plena and 
L. scutulata in some shell characters were statistically sig¬ 
nificant, no one character accurately separated individuals 
to species. Reliability of identification was enhanced by 
simultaneously using several characters through discrim¬ 
inant analysis. Table 1 lists discriminant functions for 
three analyses using the following subsets of shell char¬ 
acters: (1) length, width, presence of amber band, and 
number of whorls; (2) length, width, presence of amber 
band, number of whorls, and presence of tessellations; and 
(3) length, width, presence of amber band, number of 
whorls, and apical angle. A discriminant score was com¬ 
puted for each individual by multiplying each of the snail’s 
shell characters by the corresponding coefficient and add¬ 
ing together these products. The resulting discriminant 
score represented the number of standard deviations that 
snail was away from the mean of all snails (both species) 



DISCRIMINANT SCORE 

Figure 5 

Frequency histograms of discriminant scores for Ltttonna plena 
(n = 80) and L. scutulata ( n = 90). Discriminant score = 
-0.882L + 1.662W - 1.274B + 0.821WH + 0.473T - 
4.644 where: L = shell length in mm; W = shell width in 
mm; B = presence of amber band in aperture; YVM * number 
of whorls; T = extent of tessellated color pattern. B and T are 
scored (0) absent, (1) obscure, or (2) distinct. Distributions are 
significantly different. Kolmogorov-Smirnov two-sample test: 
DMAX = 0.7972; P c 0.01. 


'Fable 1 

Discriminant functions for three sets of shell characters: L = shell length in mm; \V = shell width in mm B = presence 
of amber band in aperture, scored 0-2; \VH = number of whorls; T = extent of tessellated color pattern scored 0 2, A = 
apical angle in degrees. Sample sizes (n) indicate the numbers of each species from open coast populations used in the 
analysis. Mean discriminant scores indicate the species’ means for each function. SI) = standard deviation. 


Littorina plena Littorina scutulata 


ysis 


Function 



n 

Mean 

SD 

n 

Mean 

SD 

1 

D = —0.732L 4- 1.599W - 

0.983 B + 1.306W11 

- 6.952 


256 

-1.128 

0.940 

215 

1.343 

1.067 

2 

D = -0.882L + 1.662W - 

1.274B + 0.821 W11 

4- 0.473T 

- 4.644 

80 

- 1.459 

0.920 

90 

1.297 

1.066 

3 

D - -1.369L + 2.514W - 

0.812B + 1.514WH 

- 0.064A 

- 3.798 

156 

— 1.119 

0.923 

125 

1.396 

1.089 
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Table 2 

Standardized discriminant function coefficients for each 
analysis of Littorina plena and L. scutulata. See text for 
explanation. 


Coefficients for each shell character 

Anal¬ 

ysis 

Shell 

length 

Shell 

width 

Amber 

band 

Whorls 

Tessel¬ 

lations 

Apical 

angle 

1 

-1.003 

1.409 

-0615 

0.629 

— 

— 

2 

-1.036 

1.299 

-0.704 

0.398 

0.331 

— 

3 

-1.991 

2.338 

-0.530 

0.702 

— 

-0.351 


on the given discriminant function. Table 1 indicates the 
number of snails used in the analyses and the mean dis¬ 
criminant scores for L. plena and L. scutulata for each of 
the discriminant functions. The extent to which the two 
species were separated by discriminant analysis is shown 
in Figure 5 (for the function that utilized shell length, 
width, presence of amber band, number of whorls, and 
degree of tessellation as discriminating variables). 

Standardized (weighting) coefficients for each discrim¬ 
inant function are given in Table 2. The absolute value 
of a coefficient indicates the relative contribution of the 
associated shell character to that particular function. (The 
sign of the coefficient simply denotes whether the char¬ 
acter is making a positive or negative contribution.) In all 
three analyses, shell width and shell length had the great¬ 
est importance in separating Littorina plena from L. scu¬ 
tulata. Apical angle and degree of tessellation were least 
important. 

Individuals were assigned to species through the use of 
classification functions. Classification functions are usu¬ 
ally employed to classify new individuals with unknown 
identities, but were used in this analysis to test the ade¬ 
quacy of the discriminant functions. As with the discrim¬ 
inant functions, classification scores were obtained by 
multiplying character values by the appropriate coeffi¬ 


cients and summing these products. For each analysis, 
there was a separate equation for each species; the re¬ 
sulting classification scores quantified the likelihood that 
an individual belonged to the corresponding species, and 
the individual was then assigned to the species for which 
it had the greatest likelihood of membership (the highest 
classification score). Classification functions and the per¬ 
cent of correct classifications ( i.e ., those confirming the 
reproductive characters) are given in Table 3 for each 
species for the three statistical analyses. (Additional spec¬ 
imens from exposed coast populations, not included in the 
initial analyses, were added to the samples of snails clas¬ 
sified.) The best classification of species was obtained when 
length, width, presence of an amber band, number of 
whorls, and presence of tessellations were employed as 
discriminating characters; better than 95% of the L. plena 
and 89% of the L. scutulata were correctly identified. 

Suitability of Derived Discriminant Functions 

Discriminant and classification functions derived for 
exposed coast populations of Littorina plena and L. scu¬ 
tulata (using length, width, presence of an amber band, 
number of whorls, and presence of tessellations as dis¬ 
criminating characters) correctly identified high propor¬ 
tions of snails collected from the same location three years 
later. The average discriminant score (mean ± SD) for 
L. plena collected at the later date was —1.492 ± 0.854 
and 97% of the individuals (n = 34) were properly clas¬ 
sified. The average discriminant score for L. scutulata was 
0.845 ± 1.092 and 82% of the individuals (n = 33) were 
properly classified. 

The same functions successfully classified 86% of the 
Littorina scutulata collected from protected harbor habitats 
(average discriminant score = 1.063 ± 0.872; n — 22). 
However, only 47% of the L. plena from protected harbor 
habitats were correctly classified (average discriminant 
score = —0.063 ± 0.971; n =57). Individuals of both 
species tended to have more whorls and a higher degree 
of tessellation when collected from the harbor compared 


Table 3 


Classification functions for each discriminant analysis. Also shown are the percentages of correct identifications for each 
species and for each analysis. (Sample sizes are given in parentheses and include additional open coast snails to those 
used in the initial analysis.) C p = classification score for Littorina plena ; C s = classification score for L. scutulata. Other 

abbreviations the same as for Table 1. 


Anal¬ 


Percent correctly identified 

ysis 

Classification functions 

L. plena 

L. scutulata 

Total 

1 

C p = — 11.76L + 24.63W + 5.07B + 18.13WH - 53.86 

C, = -13.57L + 28.58W + 2.64B + 21.36WH - 71.31 

89.82 

(275) 

88.14 

(253) 

89.02 

(528) 

2 

C p = -8.35L + 23.09W + 8.43B + 18.20WH - 0.33T - 67.71 

C, = -10.78L + 27.67W + 4.92B + 20.46WH + 0.97T - 80.29 

95.65 

(92) 

89.34 

(122) 

92.06 

(214) 

3 

C p = 38.40L - 49.21 W + 5.49B + 9.00WH + 4.54A - 174.95 

C, = 34.94L - 42.89W + 3.45B + 12.81WH + 4.38A - 184.84 

92.64 

(163) 

85.50 

(131) 

89.46 

(294) 
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to snails found on the exposed coast. These morphological 
changes in the harbor habitat increased the likelihood that 
any particular snail would be classified as L. scutulata by 
discriminant analysis, and thus significantly increased the 
probability of mistakenly classifying harbor specimens of 
L. plena (adjusted chi-square — 43.87; d.f. = 1; P < 0.001). 


DISCUSSION 

For gastropods that show as much intraspecific variability 
in shell morphology as Littonna (Whipple, 1965; Bor- 
KOWSKI, 1975; Raffaelli, 1979), morphological discrim¬ 
ination of sibling species becomes a difficult task. Mor¬ 
phological features of the shell are usually the most 
convenient, rapid, and inexpensive means for classifying 
individuals from different species, being applicable to liv¬ 
ing snails and dried shells alike. In the case of sibling 
species, however, specific diagnostic characters may be 
lacking. It may be necessary to rely upon statistical tech¬ 
niques that incorporate several shell characters in order 
accurately to identify species by morphology. 

Adult Littorina plena and L. scutulata from the exposed 
coast near Bodega Bay differ significantly in several shell 
characters, supporting previous conclusions (Murray, 
1979, 1982; Mastro et ai, 1982) that the two species are 
indeed separate taxonomic entities. Individual L. scutulata 
reach a larger size, and usually possess more whorls and 
tessellations than L. plena. Yet, classification of these in¬ 
dividuals based on shell morphology alone requires mea¬ 
surements of many morphological variables and a statis¬ 
tical procedure for evaluating these variables. Using four 
or five variables, the discriminant analyses performed in 
this study correctly assign snails to species in approxi¬ 
mately 90% of the cases. Discriminant scores further mea¬ 
sure the reliability of identification on an individual by 
individual basis; greater confidence can be placed in in¬ 
dividual identifications that have extreme discriminant 
scores. 

However, the plasticity of shell morphology in single 
species of Littorina can result in distinct morphological 
alterations along environmental and geographic gradients 
(Colman, 1932; Struhsaker, 1968; Vermeij, 1973; 
Newkirk & Doyle, 1975; Hughes, 1979). In fact, the 
shells of both L. plena and L. scutulata from Bodega Bay 
reach larger sizes, possess more whorls and tessellations, 
and have taller spires as populations occupy more shel¬ 
tered habitats. As a result, L. plena from sheltered habitats 
begin to look like L. scutulata from exposed habitats, while 
L. scutulata from sheltered habitats look like more extreme 
forms of L. scutulata from exposed habitats. Such variation 
in shell features with changes in habitat reduces the utility 
of the specific discriminant functions derived in this study. 
Clearly, systematic biases in classification occur when dis¬ 
criminant functions for exposed coast populations are ap¬ 
plied without adjustment to individuals from nearby shel¬ 
tered bays. Specimens of one species from one habitat 


overlap in morphology with specimens of the second species 
from a different habitat. 

Taxonomic separation of species based on morphology 
usually assumes that differences between species are greater 
than differences within a species. Colman (1932) states 
that two morphs are not separate species unless it is shown, 
u after the examination of sufficient numbers collected over 
a wide area, that there is not a series of overlapping in¬ 
tergrades between the two differing forms.” However, this 
stipulation for species status is conservative and makes no 
allowance for the existence of sibling species. Morpholog¬ 
ical overlap increases between the sibling species Littorina 
plena and L. scutulata when specimens are considered from 
a variety of habitats. 

Morphological characterizations of sibling species of 
Littorina are likely to lack discriminating power if consid¬ 
ered for the full geographic range of the species. Murray 
(1982) presents a discriminant analysis for L. plena and 
L. scutulata morphologies that is based on small samples 
combined from several sites along the western coast of the 
United States. Such an analysis may obscure important 
interhabitat shell variation if a wide variety of environ¬ 
ments is considered, or may seriously bias the morpholog¬ 
ical descriptions if samples from particular habitats are 
unduly represented. Snails from the exposed rocky shores 
near Bodega Bay differ significantly in morphology and 
are less accurately classified to species when methods de¬ 
veloped by Murray (1982) are utilized. The discrepan¬ 
cies between Murray’s (1982) results and the results of 
this study are likely due, at least in part, to the general 
phenomenon of variation between populations in different 
localities. In the case of sibling species, variability in mor¬ 
phological characters may be nearly as great within species 
as between species when individuals are examined over a 
wide range of habitats. 

Statistical techniques such as discriminant analysis can 
be powerful methods for distinguishing between morpho¬ 
logically similar species when other means are available 
initially to verify taxonomic status. Discriminant func¬ 
tions can be derived for particular habitats in cases where 
interhabitat variation in morphology is high, although 
other possible sources of variation (sex, season, parasit¬ 
ism) may still create difficulties in some circumstances. 
Ideally, specific techniques should be developed for spe¬ 
cific applications. Discriminant analyses performed for 
certain local populations are not suitable for species as a 
whole, nor are analyses for species as a whole likely to be 
useful for specific local populations. Discriminant func¬ 
tions are potentially important tools for discriminating 
between specimens of Littorina plena and L. scutulata and 
other sibling species of gastropods but they should be used 
with caution and frequent re-evaluation. 
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